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Abstract 

We evaluate the cross sections for the dissociation of J/ip and tj/ by tt and 
p at low collision energies, using the quark-interchange model of Barnes and 
Swanson. The dissociation cross section for J/tfi by tt is found to be relatively 
small with a maximum of about 1 mb and a kinetic energy threshold of 0.65 
GeV. The pion-induced ip' dissociation cross section is found to be much 
larger, with a maximum of about 5 mb and a lower threshold. Dissociation 
cross sections for J /ip and ip' by p mesons are also evaluated and are found 
to be large near threshold. 
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The suggestion by Matsui and Satz |l| that J /if) production might be suppressed in a 
quark-gluon plasma has led to many experimental and theoretical studies of J /if) production 
in high-energy heavy-ion collisions. The experimental observation by NA50 |2||| of anoma- 
lous J /if) suppression in Pb+Pb collisions in particular has been studied by many authors 

The evolution of a J /if) or if' produced in a heavy- ion collision depends sensitively on 
charmonium dissociation cross sections, which arise from processes such as low-energy inelas- 
tic scattering of the cc state on tt or p into open-charm final states. Small J /if) dissociation 
cross sections by tt or p may favor an interpretation of the Pb+Pb data in terms of the 
production of a new phase of matter, possibly the quark-gluon plasma. In contrast, a large 
p+J/ip dissociation cross section might imply that p+J/if inelastic scattering may be an 
important part of the Pb+Pb anomaly because the density of p mesons increases approxi- 
mately quadratically as the density of pions increases. In view of the importance of these 
dissociation cross sections for the interpretation of heavy-ion collisions, they should be eval- 
uated and incorporated in Monte Carlo simulations before any final conclusions are reached 
regarding the underlying physics. 

The dissociation of J /if by light mesons has been considered previously by several groups 



12Hl6| • Unfortunately, the numerical cross sections quoted in these references span a consid- 



erable range, due largely to different assumptions regarding the dominant scattering mech- 



anism. 



Kharzeev, Satz, and collaborators ||T2]JT3[| used the parton model and perturbative QCD 
"short-distance" approach of of Bhanot and Peskin fL7|J18|l , and found remarkably small low- 
energy cross sections for J/if> on light hadrons. For example, their J/if+N cross section at 
y/s = 5 GeV is only about 0.25 pb ||12|| . A finite-mass correction increases this cross section 



by about a factor of two However, in high-energy heavy-ion reactions, the collisions 

between the produced tt and p with J /if) and if' occur at low energies (of the order of a 
few hundred MeV to about 1 GeV relative kinetic energies). The applicability of the parton 
model and pQCD for reactions at this low energy region is open to question. 



Matinyan and Miiller |14], Haglin p5|J19|] , and Lin |20| recently reported results for these 



dissociation cross sections in meson exchange models. They use effective meson Lagrangians 
and assume t-channel D and D* meson exchange, which leads to numerical results for n+J/ip 
and p+J/ifj dissociation cross sections. Matinyan and Miiller found that these cross sections 
are rather small; both are ~ 0.2-0.3 mb at a/s = 4 GeV. Including form factors (arbitrarily 
chosen to be Gaussian with a width set to 1.5 GeV) would reduce the cross section by an 
order of magnitude. Haglin obtained a very different result, with much larger cross sections, 
by treating the D* and D* mesons as non-Abelian gauge bosons in a minimally-coupled 
Yang-Mills meson Lagrangian. Form factors were introduced in later calculations and the 
mb-scale cross sections are sensitive to the choices of the form factors pO|JT9[| . Of course 
the use of a Yang-Mills Lagrangian for charmed mesons has no a priori justification, so the 
crucial initial assumption made in these references would require independent confirmation. 
In any case, the assumption of t-channel exchange of a heavy meson such as a D or D* 
between tt and J/tp with point-like couplings is difficult to justify because the range of these 
exchanges (1/M « 0.1 fm) is much smaller than the physical sizes of the initial 7r and J /if) 
mesons. 

Charmonium dissociation processes can presumably be described in terms of the funda- 
mental quark and gluon interactions, but are of greatest phenomenological interest at energy 
scales in the resonance region. For this reason, we advocate the use of the known quark-gluon 
forces to specify the underlying scattering amplitude, which must then be convolved with 
explicit nonrelativistic quark model hadron wavefunctions for the initial and final mesons. 

Martins, Blaschke, and Quack previously reported dissociation cross section cal- 



culations using essentially the approach we describe. The short-distance interaction used 
by these authors in particular is quite similar to the form we employ. For the confining 
interaction, however, they used a simplified color-independent Gaussian potential between 
quark- ant iquark pairs only, rather than the now well-established linear X(i) ■ X(j) form. 
They found a rather large 7r+J/ip dissociation cross section which reached a maximum of 
about 7 mb at the kinetic energy in the center-of-mass system Eke of about 0.85 GeV. 



Although our approach is very similar to that of Martins et al, our final numerical results 
differ significantly, due mainly to the modeling of the confining interaction. 

In this paper we use the approach discussed above to evaluate the dissociation cross 
sections of J/ip by ir and p, and compare our results to other theoretical cross sections 
reported in the literature |I2Hni. We also calculate cross sections for the dissociation of i/j' 



by 7r and p, which have not been evaluated elsewhere. 

We employ the Barnes-Swanson quark-interchange model to determine these dis- 

sociation amplitudes. This approach uses the nonrelativistic quark potential model and its 
interquark Hamiltonian to describe hadron-hadron interactions and therefore implicitly in- 
corporates the successes of the quark model in describing the hadron spectrum and many 
static properties of hadrons. The model parameters are fixed by fits to the meson spectrum, 
so there is little additional freedom in determining scattering amplitudes and cross sections. 
One proceeds by calculating the scattering amplitude for a given process at Born order in the 
interquark Hamiltonian. In the case of meson-meson scattering, this scattering amplitude is 
given by the sum of the four quark line diagrams shown in Fig. 1. These are evaluated as 
overlap integrals of quark model wavefunctions, using the "Feynman rules" given in App. 
C of Ref . |2]J . This method has previously been applied successfully to the closely related 



no-annihilation scattering channels I = 2 hit ||21|| , 1 = 3/2 ||23|| , I = 0, 1 S-wave KN 



scattering [24| and the short-range repulsive NN interaction (25 



Following Ref. [p2j] , the interaction between each pair of constituents i and j is taken to 

be 

A(j) 

Hij ~ • ~ {^ / colorCoulomb(^?j) ^linear ('"ij) ^spin— spin yij ) ~t~ ^con} 

= m AO) | a. _ 3 _ Srcc^ f <Q ^ \ 

2 2 [nj 4 3 3m imj j \ti 3 / 2 J J V; 

This Hamiltonian is derived in the Coulomb gauge, which is the most convenient gauge for 

bound states and low-energy phenomena. The model parameter a s is the strong coupling 

constant, b is the string tension, rrii and rrij are the interacting quark or antiquark masses, 

and a is a range parameter in the Gaussian-smeared spin-spin hyperfine interaction. A 



constant shift V con is also included in the interaction. For antiquarks the generator A/2 is 
replaced by — A T /2. 




Tl T2 

Fig. 1. Born-order quark line diagrams. 

For example, a specific channel is A = J/if>, B — 7r + , 

C = D + and D = D*. 

The model parameters we employed were a s = 0.58, b = 0.18 GeV 2 , o = 0.897 GeV, 
m u = md = 0.345 GeV, m c = 1.931 GeV and V con = —0.612 GeV. This set of parameters 
gives masses within 0.08 GeV of experiment for the n, p, D(1869), L>*(2010), J/ip, and 
if)' mesons and also provides a very good description of the 1 = 2 S-wave tttt phase shift. 

An alternative set of parameters, found by fitting a large set of experimental masses, is 

i 

a s = 0.594, b = 0.162 GeV 2 , a = 0.897 GeV, m u = m d = 0.335 GeV and m c = 1.6 GeV. 
This second set, with a flavor-dependent V con , was used to test the sensitivity of our results 
to parameter variations. 

Before proceeding to our results, we note that the well-known "post-prior ambiguity" 
arises in calculations of bound state scattering amplitudes involving rearrangement collisions 
||26|| . Since the Hamiltonian which describes the scattering process AB — > CD can be 
separated into free and interaction parts in two ways, H = Hf + Hf + V AB or + 
Hp +Vcd, there is an ambiguity in the choice of Vab or Vcd as the interaction Hamiltonian. 
The first version is known as the "prior" form and leads to the scattering diagrams of Fig. 
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1, in which the interactions occur before quark interchange. The second choice is the "post" 
form, which leads to diagrams in which the interactions occur after quark interchange. 
One may show that the post and prior expressions for the scattering amplitude are equal, 
provided that exact eigenf unctions of the free Hamiltonians are used for the asymptotic states 
(The relevance of this to time reversal invariance is demonstrated numerically in Ref. 
22fl.) In our calculations we employ numerically determined Hamiltonian eigenfunctions 



for each of the external meson states considered; in the nonrelativistic case this would 
suffice to eliminate the post-prior discrepancy. In the processes considered here we have 
used relativistic kinematics and phase space, but use Galilean boosts for the states, as 
appropriate for a nonrelativistic quark model calculation. In consequence we find that the 
post and prior scattering amplitudes differ slightly. (We note in passing that one could carry 
out a relativised version of this calculation, although the full relativistic boosts would induce 
small additional effects due to Wigner rotations and creation of quarks and gluons.) In this 
paper we use the mean of the post and prior results as our theoretical cross section, and the 
estimated errors due to the post-prior discrepancy and parameter variations are indicated 
by bands in the figures. 

The cross sections we obtain for the dissociation of J / ip and ip' by 7r are shown in Fig. 2 
as a function of the kinetic energy in the center of mass system, Eke = — Ma — Mb, 
where Ma and Mb are the rest masses of the colliding particles in the initial channel. The 
lowest-lying allowed final states are D*D, DD*, and D*D* , and the total dissociation cross 
section is taken to be the sum of these three channel cross sections; this is shown as a 
solid line in the figure. ( The reactions n+J/ip — > DD and ir+ip' — > DD are AS" ^ 
transitions allowed in QCD but have zero transition matrix elements in our Hamiltonian 
(1). These transition amplitudes would be nonzero for example if we included spin-orbit 
terms in (1). The relatively weak process n+J/ip — > DD has been considered in a Dyson- 
Schwinger formalism by Blaschke et al. ||27|| , who find a maximum cross section of about 
0.1 mb near threshold. Note that the S-wave to S-wave transition is absolutely forbidden, so 
although ir+ip' — > DD is actually exothermic, it does not lead to a divergent cross section 
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at threshold.) 

The n+J/ip dissociation process is endothermic and requires an initial kinetic energy of 
0.65 GeV. The cross section shows a rapid rise above threshold (as expected for an S-wave 
process) and has a broad maximum of about 1 mb not far above threshold (Fig. 2a). This 
is somewhat smaller than the ~ 7 mb estimated by Martins et ah, which we discuss below. 



Jt + /A|/ ->D+D orD +D 

K + J/y->D* + D (a) 




E m -M % -M { ^ r} (GeV) 



Fig. 2. Cross sections for pion-induced dissociation of J/ip (Fig. 2a) and 
~ ip' (Fig. 26). The solid curves give the total dissociation cross section. Es- 
timated systematic errors due to parameter uncertainties and the post-prior 
discrepancy are shown as bands. 

The cross section for dissociation of the ip' by n is rather larger in part because this 
reaction is only weakly endothermic; the initial n+ip' kinetic energy in n+ip' — > DD* and 
DD* is only about 0.05 GeV at threshold. The total cross section reaches a maximum of 
about 6.2(0.8) mb at the kinetic energy of about 0.1 GeV and has a secondary maximum 
of 4.6(1.8) mb at the kinetic energy of about 0.22 GeV due to the opening of the D*D* 
channel. Notice that the ratio of the peak values of the ir+ip' and ir+J/ip cross sections is 
roughly 6; this should be contrasted with the prediction of ~ 5000 given in Ref. [17]. The 
minimum in the cross section near the kinetic energy of 0.4 GeV is due to the complete 
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destructive interference between transfer (Tl and T2) and capture (CI and C2) diagrams. 

We next calculate the p+J/tp and p+ip' dissociation cross sections. The allowed low- 
lying final states are DD, DD* and D*D (S tot = 0,1), and D*D* (S to t = 0,1,2). These 
cross sections are shown in Fig. 3. Since the reaction p + J/ip — > DD is exothermic, this 
cross section diverges as l/\v p j/^\ near threshold. For other channels the thresholds occur at 
higher energies, so those subprocesses are endothermic. The total dissociation cross section 
is shown as a solid line in Fig. 3. It is numerically about 11(3) mb at a kinetic energy of 0.1 
GeV, decreasing to 6(2) mb at a kinetic energy of 0.2 GeV. 



In the case of p+ip f dissociation, all the channels we consider are exothermic, so the 
low-energy divergence is quite pronounced. Our numerical results for these cross sections 
are shown in Fig. 4. The total cross section decreases from 15(2) mb to 6(2) mb as the 




30.0 



E CM -M p -M J/v (GeV) 



Fig. 3. Total and individual channel cross sections J/ip dissociation by p. 
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kinetic energy increases from 0.1 to 0.2 GeV. 

30.0 p| — ■ . ■ ■ 1 — 



25.0 




E CM -M p -M v ,(GeV) 

Fig. 4. Total and individual channel cross sections for ip' dissociation by p. 
We previously noted that our n+J/ip cross section is considerably smaller than the esti- 
mate of Ref. ||16|| , although we use a similar approach. There are several differences between 
the two approaches which lead to this discrepancy. Martins et al. assumed that the confin- 
ing interaction is an attractive Gaussian potential which acts only between quark- ant iquark 
pairs. The neglect of the quark-quark and antiquark-antiquark confining interaction amounts 
to discarding the transfer diagrams (Tl and T2) for the confining potential. Since we find 
that the transfer and capture diagram confinement contributions are similar in magnitude 
but opposite in sign (due to color factors), Martins et al. did not include an important 
destructive interference. Their use of a Gaussian rather than a linear potential will obvi- 
ously lead to quantitatively different cross sections. Furthermore, the cross section values 
are quite sensitive to the parameters used when Tl and T2 are not included. All these fac- 
tors contribute to the differences between our results and those of Martin et al. for Tt+J/ip 
collisions. 

The destructive interference between transfer and capture diagrams with spin- 
independent forces (color Coulomb and confinement) has been noted previously (see, for 
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example, Refs. pT| , f2"2"f and references cited in [p5fl). This interference explains the well- 
known spin-spin hyperfine dominance in light hadron scattering in channels such as / = 2 
7T7T and NN. In the presence of heavy quarks, however, most hyperfine interaction diagrams 
are suppressed by the large charm quark mass; this is the reason we included the color 
Coulomb and confining interactions in the present analysis. It is interesting in retrospect to 
examine the different channels and determine which of the various interactions dominates 
the amplitude. We find that the hyperfine interaction still dominates the n+J/ip dissoci- 
ation amplitude, whereas the linear confining interaction dominates n+ip', p+ip, and (for 
the kinetic energy less than 0.3 GeV) p+ip'- Above 0.3 GeV, the color Coulomb interaction 
dominates p+ijj' scattering. (We caution the reader that this decomposition depends on the 
choice of the post or prior form for the T-matrix; the results we quote are for the prior form, 
involving the diagrams of Fig. 1.) 

There is no direct experimental measurement of these cross sections to which we can 
compare our results. We found a small n+J/ip cross section which starts at a high threshold 
and a large n+ifj' cross section which starts at a low threshold. If these cross sections 
are folded with a distribution of pions with an average kinetic energy of about 200 MeV, 
we would obtain cr e g-ective( 7r + J/i>) « cr cflfectivc( 7r + V'O) which is consistent with earlier 
observation in a model of J/ip and ip' suppression in O+A, and S+U collisions @]. Hopefully, 
future Monte Carlo simulations of the dynamics of charmonium in heavy-ion collisions will 
lead to a more direct comparison. The large p+J/ip and p+ip' cross sections we have 
found imply that both the J/ip and ip' will be quickly dissociated if there is a significant 
p meson population in the medium. Since our results on the divergence of p+J/ip and 
p+ip' dissociation cross sections at threshold follow directly from simple kinematics, these 
results must be qualitatively correct. The normalization of these cross sections, however, 
required detailed calculation and should also be compared to experiment if possible. Since 
dissociations of J/ip by tc and p populate different states (for example n+J/ip doe not 
lead to DD but p+J/ip does), it may be possible to separate these processes and their 
associated cross sections by studying the relative production of DD, D*D+h.c. and D*D* 
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if the expected open charm background can be subtracted. 

In the future it may be useful to carry out detailed simulations of J/ip absorption in 
heavy-ion collisions using the cross sections obtained here to test the accuracy of our results. 
If our cross sections do prove to be reasonably accurate, it will clearly be useful to incorporate 
them in simulations of J/ip suppression in Pb+Pb collisions and in other processes that use 
charmonium as a signature of the quark-gluon plasma in order to subtract the effects of J/ ip 
suppression due to its interaction with hadron matter. 
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